Optical Properties of Correlated Materials — 
or Why Intelligent Windows may look Dirty* 



OS 

o 
o 

(N 



O 

o 



> 
in 
i> 
in 

o 

o> 
o 



X 



Jan M. Tomczak 1,2 and Silke Biermann 3 

Research Institute for Computational Sciences, AIST, Tsukuba, 305-8568 Japan 
2 Japan Science and Technology Agency, CREST 
3 Centre de Physique Theorique, Ecole Poly technique, CNRS, 91128 Palaiseau Cedex, France 

Materials with strong electronic Coulomb correlations play an increasing role in modern mate- 
rials applications. "Thermochromic" systems, which exhibit thermally induced changes in their 
optical response, provide a particularly interesting case. The optical switching associated with the 
metal-insulator transition of vanadium dioxide, for example, has been proposed for use in numer- 
ous applications, ranging from anti-laser shields to "intelligent" windows, which selectively filter 
radiative heat in hot weather conditions. 

Are present-day electronic structure techniques able to describe, or - eventually even predict - 
such a kind of behavior ? How far are we from materials design using correlated oxides ? These are 
the central questions we try to address in this article. 

We review recent attempts of calculating optical properties of correlated materials within dynam- 
ical mean field theory, and summarize results for vanadium dioxide obtained within a novel scheme 
aiming at particularly simple and efficient calculations of optical transition matrix elements within 
localized basis sets. 

Finally, by optimizing the geometry of "intelligent windows" , we argue that this kind of technique 
can in principle be used to provide guidance for experiments, thus giving a rather optimistic answer 
to the above questions. 

PACS numbers: 71.27. +a, 78.20.-e, 71.30.-rh, 71.10.-w 



TOWARDS MATERIALS DESIGN USING 
CORRELATED ELECTRON SYSTEMS ? 



Nowadays, much effort is put into finding sustainable 
ways to increase the energy efficiency of man-made pro- 
cesses. Besides classical engineering, it is materials sci- 
ence that has the potential to vastly contribute to reduce 
energy spendings. As a prototypical example, we shall 
describe a proposal made by materials scientists geared 
at saving air-conditioning costsii^. 

As a function of temperature, vanadium dioxide (VO2) 
undergoes a metal-insulator transition^. The changes 
through this transition are such that, coating a window 
with a thin layer of VO2 may - under certain conditions 
(see below) - create a fenestration that is "intelligent" 
in the following sense : At low temperatures, when VO2 
is in its insulating regime, the transmission properties 
change only remotely with respect to the bare glass win- 
dow, whereas infrared radiation is filtered when the ox- 
ide layer switches to its metallic phase in hot weather 
conditions. Then, heat radiation is substantially pre- 
vented from entering the building and there is less need 
for energy-intensive cooling by air-conditioning. 

In geographical regions, where seasonal temperature 
changes are important (and this encompasses in par- 
ticular most of the industrialized world), these "ther- 
mochromic" windows are clearly superior to "static" win- 



*This work previously appeared as the "Scientific Highlight of the 
Month" No. 88, August 2008 of the * fe -network. 



dow coatings^, that filter certain wavelengths irrespective 
of external conditions and which thus would increase the 
need for heating buildings during winter. 

Even if in vanadium dioxide the changes of its proper- 
ties with temperature are particularly abrupt, VO2 is by 
no means the only material where a tiny change in exter- 
nal parameters can radically modify the physical prop- 
erties. In fact, such behavior can nearly be considered a 
hallmark of materials with strong electronic Coulomb in- 
teractions. The concerted behavior of electrons in corre- 
lated materials causes indeed quite in general an extreme 
sensitivity to external stimuli, such as temperature, pres- 
sure or external fields. Heating insulating SmNi03 be- 
yond 400 K or applying a pressure of just a few kbar to 
the Mott insulator (Vi-xCr^Os (x=0.01)&, for exam- 
ple, makes the materials undergo transitions to metallic 
states. This tuneability of even fundamental properties 
is both, a harbinger for diverse technological applications 
and a challenge for a theoretical description. 

An increasing role is nowadays played by artificial 
structures, ranging from multilayers that display the gi- 
ant magnetoresistance effect (widely used in storage de- 
vices to functional surfaces, where appropriate coat- 
ings e.g. provide a self-cleaning mechanism^. The huge 
freedom in the design, which concerns not only parame- 
ters such as the chemical composition, the doping and the 
growth conditions, but also the geometry of the device 
(e.g. the layer thicknesses), however makes the search for 
devices with specific electronic properties a tedious task. 
This leads us to the central question of the present arti- 
cle : Can modern first principles calculations help in the 
quest for promising materials and setups for particular 
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devices ? The aim can of course not be to replace exper- 
iments, but rather to provide some guidance in order to 
minimize expensive experimental surveys and prototyp- 
ings. 

The Achilles heel of electronic structure theory is the 
description of electronic many-body interactions. Indeed, 
a great majority of materials used in modern applica- 
tions fall in the class of so-called "strongly correlated ma- 
terials" where electron-electron interactions profoundly 
modify (if not invalidate) a pure band picture^. State- 
of-the-art first principles methods, such as density func- 
tional theory (DFT) 1 ^, are then no longer sufficient to 
predict the physical properties of these materials. De- 
spite these difficulties, in this article we give a quite op- 
timistic view on the above question. 

In fact, important steps to bridge the gap between 
band structure methods and many-body physics have 
been made in recent years. For materials with moder- 
ately strong correlations, Hedin's GW approximation- 
- which has seen increasingly sophisticated implementa- 
tions in the electronic structure contex t 12 i 13 ' 14 i 15 - has 
established itself as a method of choice. For strongly 
correlated materials, progress was for instance brought 
about by combining density functional theory within 
the local density approximation 1 ^ with dynamical mean 
field theory (DMFT)^ 7 .. The resulting approach, dubbed 
LDA+DMFTiSiiS (for reviews see 2 ^ 2 ! ^ ) , joins the ac- 
curate description of strong local Coulomb correlations 
in a many-body framework with the material-specific in- 
formation provided by state-of-the-art band theory. It 
has - over the last years - helped to elucidate physical 

mechanisms at work in systems such as transition met- 
als 23 i 24 i 25 1 26 j their 0^^27,28,29,30,31,32,33,34,35,36 or sul- 
phides 3 ^! 3 ^, as well as f-electron compound o 39 ' 40 i 41 ' 42 . 
The number of applications is nowadays too large to give 
a complete list in this work, but recent reviews provide 
an extensive pictur o 20 i 21 ' 22 i 43 . 

While tremendous progress has been achieved, there 
remains a chasm between what state-of-the-art electronic 
structure methods can calculate and what experimental- 
ists are actually measuring. On the one hand, the chem- 
ical complexity of many systems does simply not (yet) 
allow for being tackled by costly many-body techniques. 
On the other hand, the variety of experimental observ- 
ables that are being computed for the sake of comparison 
is rather unsatisfactory. It is therefore an important task 
to make more experimentally measurable quantities ac- 
cessible from theoretical calculations. 

Within dynamical mean field theory, emphasis is com- 
monly put on spectral properties, and the evaluation of 
observables other than spectral functions is a rather new 
advancement in the realistic context. Yet, it is rather the 
response behavior of correlated materials that is promis- 
ing for applications. 

In this article, we review some recent attempts of cal- 
culating optical properties of correlated materials, and 
explore the implications for technological applications 
on the specific example of V02-based intelligent win- 



dows. For further reading and some of the original work 

seQ 44^16 



II. OPTICAL SPECTROSCOPY 

A. General formulation and some physical 
implications 

Numerous experimental techniques have been devised 
for and applied to the study of correlated materials of 
ever growing complexity. Optical spectroscopy, which 
is the subject of this work, is, in a way, the most nat- 
ural among them : Optical detectors are sampling the 
response to incident light, as do our eyes, albeit access- 
ing frequencies, and thus phenomena, that are beyond 
our vision. The technique is particularly suited to track 
the evolution of a system under changes of external pa- 
rameters like temperature or pressure. This is owing to 
a generally high precision, and the fact that, contrary 
to e.g. photoemission spectroscopy or x-ray experiments, 
results are obtained in absolute values. Especially, the 
existence of sum-rules (see e.gi 47 i 48 ) allows for a quanti- 
tative assessment of transfers of spectral weight upon the 
progression of the system properties. Moreover, while in 
photoemission the electron escape depth and thus sur- 
face effects are often an issue, the larger skin penetration 
depth assures that optical spectroscopy is a true bulk 
probe. One might add that the transition matrix ele- 
ments are also better understood in optics (this is an 
important part in this work) than in photoemission (see 
e.gji 9 -). On the other hand, response functions are two- 
particle quantities that are less obvious in their inter- 
pretation than a one-particle spectrum. An important 
simplification is achieved when neglecting vertex correc- 
tions. In the framework of linear response theory, the op- 
tical conductivity can then be expressed as (for reviews 
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Here, A^(lu) is the momentum-resolved many-body 
spectral function, and different optical transitions are 
weighted by the Fermi velocities v^ a = — (kL'\T a \kL), 
matrix elements of the momentum operator V. Both, 
spectral functions and velocities are matrices in orbital 
space L, which we will specify later on. The Fermi func- 
tions i(cu) select the range of occupied and empty en- 
ergies, respectively, V is the unit-cell volume, a,/3 de- 
note cartesian coordinates, and Re a a0 is the response in 
a— direction for a light polarization E along j3. 

To get an idea about the physical content of this for- 
mula, we first consider some archetypical cases. At zero 
temperature, and for a system that is well described by 
its band structure as given by density-functional theory 
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based methods, the (Kohn-Sham) spectrum is the defin- 
ing quantity, since vertex correction are absent^!. More- 
over, the spectral functions A(k, lu) of the system become 
Dirac distributions and the trace in Eq. |T|) reads2& 

5{uj + lu + fi-e k )v k <H^ +A*-£k)«k ( 2 ) 

m,n 

While at finite frequencies, lu > 0, only inter-band transi- 
tions, m 7^ n, can give a contribution, we see that at zero 
frequency, lu = 0, and provided a band e£ is crossing the 
Fermi level, the response is a delta function deriving from 
intra-band transitions^. The latter is just the result of 
the fact that without interactions (electron-electron cor- 
relations or a coupling to a bosonic mode) or disorder, 
the lattice momentum k is a constant of motion and the 
current thus does not decay. 

In the (effective) non-interacting case the response of 
a metal to an electric field is thus infinite. It was P. 
Drudc in 1900 who derived an expression that takes into 
account the finite lifetime of the electron excitations, by 
introducing a relaxation time r for the charge current 
j(t) = j(0)e _t / T . His expression for the conductivity can 
be recovered from our linear response result, Eq. ([1]), by 

f-2 jr 2 

assuming free particles, e k = and a constant self- 

energy £ — — The resulting optical conductivity is 
then given by 



with n being the average charge carrier density. Hence, 
the response at zero frequency is finite and non-zero at 
finite energies. 

While the physical picture about the origin of the re- 
laxation time had to be revised with the advent of quan- 
tum mechanics and Bloch's theory of electronic states in 
periodic potentials, Drude's theory still accounts well for 
the response of simple metals. However, it becomes insuf- 
ficient, even in the one-band case, in the presence of sub- 
stantial correlation effects. Indeed, the influence of elec- 
tronic interactions is beyond a mere broadening of bands. 
For the case of the one-band Hubbard model within 
DMFT, the paramount characteristic is the appearance 
of Hubbard satellites in the spectral function. Therewith, 
not only transitions within a broadened quasi-particle 
peak are possible (a la Drude), but also transitions from 
and to these Hubbard bands aris o 51 ' 52 . Thus, in such a 
metal two additional contributions occur, stemming from 
transitions between the quasi-particle peak and the indi- 
vidual Hubbard bands. This feature is often referred to 
as the mid-infrared peak, due to its location in energy in 
some compounds. At higher energy, transitions between 
the two Hubbard bands appear. In the Mott insulating 
phase, only the latter survive. Therewith the complex- 
ity of optical spectra is considerably enhance d 51 i 52 with 
respect to the independent particle picture. 

Within DMFT, calculations of the optical conductivity 
were first performed by Pruschke and Jarrell et al. 53 i 54 



for the case of the Hubbard model. Rozenbcrg et al. 51 i 52 
studied the phenomenology of the different optical re- 
sponses of the Hubbard model throughout its phase dia- 
gram, discussing the above described phenomenology and 
comparing it to experiments on V2O3. With the advent 
of LDA+DMFT, optical conductivity calculations gained 
in realism: Bliimer et al. 55 i 56 and later on Pavarini et al. 
■21 and Baldassarre et al. ^ used the LDA+DMFT spec- 
tral functions for the calculation of titanate and vanadate 
optical spectra. A more general approach was developed 
by Palsson^ for the study of thermo-electricity. Our 
wor k 36 i 44 i 45 i 46 goes along the lines of these approaches. 
We will however employ a formulation using the full va- 
lence Hamiltonian therewith allowing for the general case 
including intcrband transitions and we extend the inter- 
vening Fermi velocities to multi-atomic unit cells, which 
becomes crucial in calculations for realistic compounds. 

Alternative techniques were proposed by Perlov et al. 
60 ' 61 and by Oudovenko et al. — . The former explic- 
itly calculated the matrix elements, albeit using a differ- 
ent basis representation than in the DMFT part, while 
the latter diagonalized the interacting system, which al- 
lows for the analytical performing of some occurring in- 
tegrals due to the "non-interacting" form of the Green's 
function. Owing to the frequency-dependence of the self- 
energy, however, the diagonalization has to be performed 
for each momentum and frequency separately so that 
the procedure may become numerically expensive. This 
technique has in particular been applied to elucidate the 
a — 7 transition in Ce^ 3 -, and, more recently, to the heavy 
fermion compound CelrlngH. 

The calculation of accurate absolute values of the op- 
tical response has however turned out to be a challenge 
which mainly stems from the way the Fermi velocities 
are treated. In view of predictive materials design ac- 
curate absolute values in the optical conductivity are a 
condition sine qua non. 

B. A scheme for the optical conductivity using 
localized basis sets 

Correlation effects enter the calculation of the opti- 
cal conductivity, Eq. |T]), only via the spectral functions 
Ak(w), while the Fermi velocities 

vt'^ = -(kL'\V a \kL) (4) 
m 

are determined by the one-particle part of the system 1 ^. 

While the computation of the latter is straightforward 
in a plane-wave setup, the application of many-body 
techniques such as LDA+DMFT, necessitates the use of 
localized basis sets (see e.g^). It is thus convenient to 
express also the Fermi velocities in terms of this basis. 
Therewith, however, the evaluation of matrix elements 
of the momentum operator, Eq. ((J), becomes rather te- 
dious. The aim of the current approach is to employ a 
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controlled approximation to the full dipole matrix ele- 
ments in a Wannier like basis set that is accurate enough 
to yield absolute values of the conductivity to allow for 
a quantitative comparison with experiment. 

We choose our orbital space by specifying L = 
(n, l,m, 7), with the usual quantum numbers (n,l,m), 
while 7 denotes the atoms in the unit cell : |ki) then 
is the Fourier transform of the Wannier function x RL (r) 
localized at atom 7 in the unit cell R. Extending the 
well-known Peierls substitution approach for lattice mod- 
els (see the review^ 7 -) to the realistic case of multi-atomic 
unit cells, we find that 
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Here, pl denotes the position of an individual atom 
within the unit-cell. The term in brackets, which is used 
in the actual calculations, is in the following referred to 
as the "generalized Peierls" term : While the derivative 
term is the common Fermi velocity, the term propor- 
tional to the Hamiltonian originates from the generaliza- 
tion to realistic multi-atomic unit-cells and accounts for 
the fact that while the periodicity of the lattice is deter- 
mined by the unit-cell, the Peierls phases couple to the 
real-space positions of the individual atoms. The correc- 
tion term that recovers the full matrix clement is denoted 
T (for its explicit form se o 36 i 45 ). The latter reduces to 
purely atomic transitions, (R, 7) = (R',7'), in the limit 
of strongly localized orbitals. In other words, the accu- 
racy of the approach is controlled by the localization of 
the basis functions. This generalized Peierls approach 
has in particular been shown to yield a good approxima- 
tion for systems with localized orbitals, such as the 3d 
or 4f orbitals in transition metal or lanthanidc/actinidc 
compounds 3 ^. 



III. APPLICATION : FROM BULK VANADIUM 
DIOXIDE TO INTELLIGENT WINDOWS 

A. A brief reminder about the electronic structure 
of V0 2 

The electronic structure of VO2 and its metal-insulator 
transition has been the subject of numerous theoretical 
studies. We give here only a brief list of prior work, for 



reviews see- 



36.(57 



In particular, we here do not enter at 
all in the decade-long discussion on whether the metal- 
insulator transition in VO2 should be considered as a 
Peierls- or a Mott transition. Our point of view on this 
question is summarized in 3 - 5 -. 

While capturing structural properties surprisingly 
well^, band-structure methods do not reproduce experi- 
mental spectra: in the metal, incoherent weight at higher 
binding energies is absent, and in the insulator the gap 
is not opened in the corresponding Kohn-Sham spectra. 
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FIG. 1: Comparison of valence band photoemission spectra 
of Koethe et al. (top) and LDA+DMFT results from^ 
(bottom). Picture frorn^. 



Nevertheless, thorough LDA studie a 67 ' 71 gave useful in- 
dications and paved the road to the application of more 
sophisticated many-body approaches. 

LDA+DMFT results for the spectral properties of VO2 
agree well with experimental findings in both, the metal- 
li o 33 i 72 ' 73 and the insulating phase 3 - 3 -, as can be seen from 
Fig. [1] which compares photoemission and LDA+DMFT 
spectra. 

Since at least the insulating phase of VO2 is in fact rel- 
atively band-like 3 - 5 -, an interesting alternative approach is 
provided by the GW approximatio n 11 ' 12 . While pioneer- 
ing early workZ 4 - had to resort to a simplified scheme, re- 
cently it became possible to perform fully ah initio GW 
calculations for VO^ 1 ^^. Full GW calculations for 
the optical conductivity and the spectral function thus 
seem to come into reach, opening the way to systematic 
comparisons with DMFT results. 

Our calculation of optical propertie s] 36 ' 44 is footing on 
the LDA+DMFT electronic structure of Refi^ 3 - and our 
recent extension thereo f 34 ' 35 . Since this many-body cal- 
culation used a downfolded Hamiltonian, we use an up- 
folding scheme to include optical transition from, to and 
between higher energy orbitals (For details see^ 6 -). 



B. The bulk conductivity 

In Fig. Ufa) we display the LDA+DMFT theoretical 
optical conductivity of the high temperature phase of 
VO2 as a function of frequency and in comparison with 
several experimental dat a 68 ' 69 ' 70 , see also 7 - 7 -. As can be 
inferred from the crystal structure^ 7 -, the optical response 
depends only weakly on the polarization of the incident 
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FIG. 2: Optical conductivity of a) metallic, b) insulating VO2 for polarizations E. Theory (red) ([aab] = [0.85 0.85 0.53]), 
experimental data (i) single crystals^ (green), (ii) thin filrn^ (solid blue), (iii) poly crystalline filrn^ (dashed blue). 



light. The Drude-like metallic response (see also above) 
is caused by transitions between narrow vanadium 3d 
orbitals near the Fermi level. As a consequence, it only 
influences the low infra-red regime a crucial observation 
as we shall see in the following. The shoulder at 1.75 eV 
yet stems from intra-vanadium 3d contributions, while 
transitions involving oxygen 2p orbitals set in at 2 eV, 
and henceforth constitute the major spectral weight up 
to the highest energies of the calculation. 

In Fig. [2fb) we show the conductivity for insulating 
VO2. As was the case before, the theoretical results are 
in good agreement with the different experiments. This 
time, a slight polarization dependence is seen in both, ex- 
periment and theory, owing to the change in crystal sym- 
metry. Indeed, vanadium atoms pair up in the insulator 
to form dimers along the c-axis, leading to the forma- 
tion of bonding/anti-bonding states for 3d orbital o 10178 . 
Optical transitions between these orbitals result in an 
amplitude of the conductivity that, in the corresponding 
energy range (uj = 1.5 — 2.5 eV), is higher for a light po- 
larization parallel to the c-axis than for other directions 
(For a detailed discussion se o 36 ' 45 ). 



C. Intelligent windows 

Having established the theoretical optical response of 
bulk VO2, and thus verified that our scheme can quan- 
titatively reproduce optical properties of correlated ma- 
terials, we now investigate the possibilities of V02-based 
intelligent window coating o 2 ' 3 ' 81 i 82 . The effect to be ex- 
ploited here can already be seen in the above responses 
of the bulk : The respective conductivities of both phases 
(Fig. Uta), (b)) exhibit a close similarity in the range of 
visible light (uj = 1.7 — 3.0 eV), whereas in the infra- 
red regime (uj < 1.7 eV) a pronounced switching occurs 
across the metal-insulator transition. As a result, heat 



radiation will be let through at low external tempera- 
tures, while its transmission will be hindered above the 
transition, which reduces the need for air conditioning 
in e.g. office buildings. The insensitivity to temperature 
for visible light, in conjunction with the selectivity of 
the response to infrared radiation, is an essential fea- 
ture of an intelligent window setup. Yet, for an appli- 
cable realization, other important requirements have to 
be met. First of all, the switching of the window has to 
occur at a relevant, i.e. ambient, temperature. Also, the 
total transmittancc of V02-films needs improvement in 
the visible ranged 3 -, and the visible transmission should 
be rather independent of the wavelength such as to pro- 
vide a colorless vision. Experimentalists have addressed 
these issues and have proposed potential solutions^ 3 - : Di- 
verse dopings, M x Vi- x C>2, were proven to influence the 
transition temperature, with Tungsten (M = W) be- 
ing the most efficient : A doping of only 6% resulted 
in T c w 20°G2i. However, this causes a deterioration 
of the infrared switching. Fluorine doping, on the other 
hand, improves on the switching properties, while also 
reducing T ^i 86 . An increase in the overall visible trans- 
mittance can also be achieved without modifying the in- 
trinsic properties of the material itself, but by adding 
antireflexion coatings with for example TiO^- 

Here, we address the optical properties of window coat- 
ings from the theoretical perspective. In doing so, we 
assume that the specular response of V02-layers is suf- 
ficiently well-described by the optical properties of the 
bulk, and we use geometrical optics to deduce the prop- 



erties of layered structures, as shown e.g. in Fig. 3(a) 



Fig. 4(a) We employ a technique equivalent to the trans- 
fer matrix method (see e.gM-), which accounts for the 
multiple reflexions within the layers. The reflectivity of 
the setup is then given by a recursion formula. Defining 
for the j th layer the phase factors (Xj(uj) — exp(iuj/cnjSj), 
with 5j and hj(ui) being the layer thickness and the 
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FIG. 3: Setup of an intelligent window, (a) geometry of the setup : VO2 on SiC>2. (b) Specular reflectivity at high (dashed) 
and low temperatures (solid). Theory (red) : 60 nm VO2 on Si02. Experiments : 50 nm VO2 on SiO-2-^ ( green), 50 nm VO2 
on Pyrex glasses (blue). All theoretical data for E || [110] polarization, and a 3mm Si02 substrate. 



complex refractive index, respectively, and denoting the 
"bare" reflection and transmission coefficients of an iso- 
lated interface between semi-infinite layers j — 1, j by 
Tj-ii, ij-i.j one obtains a recursion formula for the 
"dressed" quantities r .j, tn.j that account for the re- 
sulting properties of the stack of the layers 0, 1, • • • , j : 



?a,j — Foj-i + to,j-itj-io 



1 - r j-hJ f j-ho & j-i 



(6) 



fj,o = rjj-i +t jd - 1 t j ^ij- — '- (7) 

1 _ r 3-1-,3 r 3-1-fi a j-l 

where dj = IIfc<j a k is an effective phase, tj t k = 1 — fj 7 k 
the complex transmission, and r^o the dressed reflectivity 
for light that reaches the setup from the opposite side. 



1. A VOi-coated window 

First, we consider the most simple possible setup, 
which consists of a single V02-laycr (of thickness d VQ2 ) 
on a glass substratei£2. Such a window has been exper- 
imentally investigated by Tazawa et al. — and Jin et al. 
' 9 . In Fig. [3]we show their measured reflectivity data as 
a function of wavelength, in comparison with our theo- 
retical results. 

At low temperatures, i.e. insulating VO2; the calcu- 
lated reflectivity is in quantitative agreement with the ex- 
perimental data. In the visible range (A « 400 — 700 nm), 
the reflectivity depends strongly on the wavelength. 
Therefore, the current window will filter certain wave- 
lengths more than others, resulting in an illumination of 
a certain color - an obvious drawback. Moreover, the 
reflectivity in this region is rather elevated, causing poor 
global transmission. In the infrared regime (A > 700 nm) 
- and beyond - the reflectivity decreases, and heat ra- 



diation can pass through the setup. At high tempera- 
tures, the infrared reflectivity switches to a rather ele- 
vated value, thus filtering heat radiation. The changes 
in the visible region are less pronounced, but still per- 
ceptible, and both, the degree of transparency and the 
color change through the transition. As a result, though 
presenting qualitative features of an intelligent window, 
the current setup is not yet suited for applications. 



2. Improving the window with an anti-reflexion layer 

Having established the accuracy of our approach also 
for the case of layered structures, we now investigate a 
more complicated setup. As depicted in Fig. Ufa), an ad- 
ditional (rutile type) Ti02-coating is added on the VO2- 
layer, with the objective of serving as an antireflexion 
filter—. With the thicknesses d VQ2 and d T io 2 > the geome- 
try of the current setup thus has already two parameters 
that can be used to optimize the desired optical proper- 
ties. Since, however, each variation of them requires the 
production of a new individual sample under compara- 
ble deposition conditions, along with a careful structural 
characterization in order to guarantee that differences in 
the optical behavior are genuine and not related to varia- 
tions of the sample quality, the experimental expenditure 
is tremendous. This led Jin et al. — to first estimate a 
highly transmitting setup by using tabulated refractive 
indices and to produce and measure only one such sam- 
ple. Here, we shall use our theoretical results on VO2 
to not only optimize the geometry (d T io 2 1 4o 2 ) with re- 
spect to the total visible transmittance, but we shall be 
concerned with yet another important property of the 
window that needs to be controlled : its color. 

Fig. U displays the normalized visible specular 
transmittancei2£ for our window in its high (b) and 
low (c) temperature state, as a function of both film 
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(a)sctup of a multi-layer window (b)visible transmittance for metallic VO2 (c)visible transmittance for insulating VO2 

FIG. 4: Setup of a multi-layer intelligent window, a) geometry of the setup : Ti02 on VO2 on Si02- b), c) Setup with 
antireflexion coating : Ti02 / VO2 on SiC>2. Shown is the normalized visible transmittance (see text) an the corresponding 
color of the transmitted light for b) high and c) low temperature as a function of the layer thicknesses d v02 and d T io 2 • All 
theoretical data for E || [110] polarization, and a 3mm SiC>2 substrate. 



thicknesses. On the same graph, we moreover show the 
corresponding transmission color. The evolution of the 
light interferences within the layers results in pronounced 
changes in both, the overall transmittance and the color. 
The coating of VO2 globally degrades the transparency of 
the bare glass window. An increase of the Ti02-coating, 
on the other hand, has the potential to improve on the 
total transmittance. This can be understood from the 
mechanism of commonly used quarter- wave filters. The 
wavelength-dependence of the real-part of the TiC>2 re- 
fractive index, n TiC , 2 (A) , results in an optimal quarter- 
wave thickness, (5 Ti02 (A) = A/(4n Ti02 (A)), which varies 
from blue to red light only slightly from <5 TiC , 2 (A) = 40 
to 60 nm. This and the fact that the imaginary part 
of the refractive index, fc Ti02 (A), is negligible for visi- 
ble light also explains why the color does not change 
significantly with d Ti02 . While, as for TiC>2, the varia- 
tion of the real-part of the VO2 refractive index yields 
a rather uniform ideal thickness <5 V o 2 (A), its imaginary 
part changes significantly (by a factor of 4) within the 
range of visible light. As a consequence, the color is very 
sensitive to V02-dcposition. At higher thickness d V o 2J 
however, this dependence becomes smaller and the color 
lighter, as seen in Fig. 0Jb) and (c). Our theoretical 
transmittance profiles suggest relatively thick windows 
to yield good visual properties. Indeed, at low tem- 
peratures (Fig. 0Jc)) the local maximum that gives the 
thinnest window (Fig.^b)) is located at (d T io 2 1 ^vo 2 ) ~ 
(40 nm, 85 nm) within our calculation. However, this 
setup is still in the regime of important color oscil- 
lations. Given the uncertainties in industrial deposi- 
tion techniques, it seems rather cumbersome to consis- 
tently stabilize colorless samples. From this point of 
view, a thicker V02-film would be desirable. Indeed, 
while almost preserving the overall transmittance, a col- 
orless window at low temperatures is realized in our 
calculation for (c? T io 2 , ^vo 2 ) ~ (50 nm, 220 nm), or for 
(d Ti o 2 ,rf V o 2 ) = (> 100 nm, 220 nm). In the high tem- 
perature state, Fig. BJb), the transmittance is globally 



lower than at low temperatures. Moreover, only the 
(d Ti o 2 ,rf V o 2 ) = (> 100 nm, 220 nm) setup exhibits a 
simultaneously high transmittance in both states of the 
window. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this article, we have reviewed a scheme for optical 
properties of correlated materials that is geared at fast 
and accurate calculations within localized basis sets. 

As an example, we have addressed - for the first time 
from a theoretical perspective - the fascinating applica- 
tion of the seemingly simple composition of vanadium 
dioxide in so-called intelligent windows. 

Even if, from a technical point of view, possible ways 
of improvement in view of a fully general optics scheme 
are obvious - ideally, one may want to include e.g. ver- 
tex corrections, replace the LDA+DMFT starting point 
by a GW+DMFT calculation^ or calculate at least the 
ligand-orbital energies by many-body techniques such as 
GW - our simple and efficient optics scheme has proven 
to reach quantitative agreement with experiments, at 
least for materials where excitonic effects are negligible. 
We have shown that it is not only useful to address ques- 
tions of fundamental physics, but that one can indeed 
think of exploiting present day techniques for applied 
purposes. 

One may argue, of course, that a crucial ingredient 
entering LDA+DMFT calculations is a reliable estimate 
for the local Coulomb interactions, the Hubbard U. In 
general, it is probably a fair statement to say that - 
despite a number of techniques that have by now been 
proposed for the calculation of U (constrained LDA2 9 -, 
constrained RPA2&2I, or GW+DMFT^^^) - its de- 
termination still presents one of the bottlenecks for ab 
initio materials design. The recently much discussed new 
family of iron pnictidc superconductors provides an in- 
teresting test case for calculations of U 94 i 95 i 96 , where we 
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dispose of little a priori knowledge. 

For vanadium dioxide, however, this issue seems to be 
more a conceptual than a practical one, since experimen- 
tal and theoretical estimates finally converge towards a 
common answer. 

In conclusion, we thus give an optimistic answer to 
the central question of this article : electronic structure 
techniques can - if not design - at least help guiding the 
search for functional materials and their devices, and this 
even for the particularly challenging class of correlated 
materials. 
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uses spectral functions stemming from a many-body cal- 
culation, (b) The particle-hole interaction (vertex correc- 
tions). This is not addressed in the above formalism and 
remains a challenge for future work. We note that vertex 
corrections vanish within DMFT (infinite dimensions) in 
the one-band case—. For the current multi-orbital case, 
this is however an approximation. Within GW— , vertex 



corrections are included on the RPA level. However, GW 
spectra may not be sufficient for strongly correlated sys- 
tems. 

for simplicity we assume to work in the Kohn-Sham basis, 
i.e. the spectral function A^uj) is diagonal, 
or transitions within degenerate bands, e£ = while 
ra 7^ n. 

While for lattice models this is only true for local density- 
density interactions, in the continuum formulation of a 
solid this requires the interaction only to be local and of 
two-body kind. Then the interaction part of the Hamilto- 
nian H = Ho + Hint commutes with the position operator 
TZ a and thus V a = —vm/h [7£ Q ,Ho]. For details see—. 
We stress that the bonding/anti-bonding splitting is not 
sufficient to open a gap within band-structure approaches. 
Indeed correlations enhance this splitting (after having 
transfered the charge into the a\ g band), and we refer 
to this scenario as a "many-body Peierls insulator"—. 
We assume the bare window to be made of quartz glass, 
SiC>2. All auxiliary refractive indices (other than the one 
of VO2) are taken from Handbook of Optical Constants of 
Solids by Edward D. Palik, Academic Press, 1985. 
floon™ d\S(X)T(X)/ f™°™ d\S(X), with the spectrum 
of the light source S(X) , and the transmittance T = 1 — R, 
R being the specular reflexion given by Fresnel's formulae 
(see e.g.—). We thus neglect absorption due to inhomo- 
geneities that lead to diffuse reflexion. This is justified 
for our applications to windows. Also, VO2 has a glossy 
appearance and hence a preponderant specular response. 



